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Abstract. The J-matrix method is extended to difference and g-difference operators and 
is applied to several explicit differential, difference, g-difference and second order Askey- 
Wilson type operators. The spectrum and the spectral measures are discussed in each case 
and the corresponding eigenfunction expansion is written down explicitly in most cases. In 
some cases we encounter new orthogonal polynomials with explicit three term recurrence 
relations where nothing is known about their explicit representations or orthogonality mea- 
sures. Each model we analyze is a discrete quantum mechanical model in the sense of Odake 
and Sasaki [J. Phys. A: Math. Theor. 44 (2011), 353001, 47 pages]. 
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1 Introduction 

The J-matrix method started with the pioneering works of Yamani, Heller and Reinhardt 
[19, 20, 40] in the early 1970's and has been applied by Yamani, Heller and Reinhardt to 
different physical models. Some of the recent applications of the J-matrix method to physics 
are spearheaded by Alhaidari and his research team, [2, 3, 4, 5, 7]. The J-matrix principle says 
that the spectrum of a tridigonalizable operator is the same as the tridiagonal matrix repre- 
senting it. Such a tridiagonal matrix can be split into irreducible blocks, and to each of these 
blocks there is a corresponding set of orthogonal polynomials. Moreover, the eigenfunctions of 
a tridiagonalizable operator can be expressed using these corresponding polynomials, and in the 
self-adjoint case the spectral measure is related to the orthogonality measures of the orthogo- 
nal polynomials. This general set-up is described and proved in [24] where we considered the 
Schrodinger equation with the Morse potential as an example. Our later work [23] develops 
a general scheme for tridiagonalizing differential, difference, or g-difference operators arising 
from two sets of related orthogonal polynomials. In particular, we find in [23] the Jacobi trans- 
form and its special case the Mehler-Fock transform originally introduced by Mehler to study 
electrical distributions. 

Tridiagonalization of an explicit symmetric or self-adjoint operator T, like a differential or 
(g-)difference operator on an explicit Hilbert space Ti of functions, amounts to finding an explicit 
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orthonormal basis {e„},^Q so that we can realize 

Ten = Oinen+l + fine-n + On-ien, > 0, /3„ G M, 

with respect to this basis. Since three-term operators, or Jacobi operators, can be described in 
terms of orthogonal polynomials, we can obtain information on T in terms of properties of the 
orthogonal polynomials, see [24] and references given there. Assuming that there are orthgonal 
polynomials pn (corresponding to a determinate moment problem) satisfying 



then the spectral decomposition is UT = MC/, where U : H ^ -^2(/^)i ^ Pn and M : L2{fi) — )• 
L2{p) is the multiplication operator. The link between more general differential, difference or q- 
difference operators, and Jacobi matrices can be very useful to study the spectral decomposition 
of the original operator in terms of the orthogonal polynomials and vice versa. So we can obtain 
information on one of the operators by transferring information from the other, and we show 
this in particular examples in this paper. In particular, in case information on both operators is 
known, we obtain even more explicit results, and examples of this approach can be found in [23]. 
However, it is not straightforward to find explicit tridiagonalization of a given operator, and here 
we present ways to obtain operators with a tridiagonalization. In the tridiagonalization of the 
operators in this paper it is often the case that the polynomials cannot be matched directly with 
known polynomials [21, 29], and in these cases we have given some information on the support 
of the spectral meaure of these polynomials. A treatment of the spectral theory of differential 
operators can be found in many sources, and we refer the interested reader to the excellent book 
by Edmunds and Evans [15]. The spectral theory of tridiagonal matrices and their connection 
with orthogonal polynomials and the moment problem is in [1, 14, 30, 38]. 

The contents of the paper are as follows. In Section 2 we record the formulas used in the se- 
quel. The basic references are [8, 16, 17, 21, 29, 37]. The expert reader may easily skip Section 2. 
In Sections 3 and 4 we start with an operator with known orthogonal polynomial eigenfunctions 
then multiply it by a linear function of the space variable and study the spectral properties of 
the new Hamiltonian. Section 3 treats the case of Laguerre polynomials, leading to tridiagona- 
lization involving continuous dual Hahn polynomials. It is simple enough but contains all the 
ingredients of the method. Section 4 treats the Meixner polynomials and the J-matrix method 
leads to a one parameter generalization of the continuous dual Hahn polynomials. The examples 
in Sections 3 and 4 are related to the approach of [23]. In Section 5 we introduce a different 
modification. We start with an operator T which is diagonalized by a known polynomial system. 
We then consider the Schodinger operator T + ^x. Such an equation arises for example in the 
case of a charged particle in the presence of a uniform electric field. In this case 7 = —qF, 
where F is the magnitude of the electric field and q is the electric charge, see [32, § 24], or [10]. 

The Schodinger operator T + jx is automatically tridiagonal in the orthogonal polynomial 
basis which diagonalizes or tridiagonalizes T. We study the spectral properties of T + 73; for four 
different sets of polynomials. In Section 6 we combine both generalizations of Sections 3 and 4 for 
the case of the Laguerre and Meixner polynomials. Finally, in Section 7 we study this approach 
for two families of g-orthogonal polynomials, namely the Al-Salam-Chihara polynomials and 
the (^"^-Hermite polynomials. The g~^-Hermite polynomials correspond to an indeterminate 
moment problem, so we study the corresponding q'-difference operator on the weighted L^-space 
corresponding to a A^-extremal measure. It turns out that the polynomials in the tridiagonali- 
zation is again corresponding to an indeterminate moment problem, so that the g-difference 
operator is not essentially self-adjoint on the space of polynomials. 

We end by noting that more differential, difference and g-difference operators can be studied 
using the J-matrix method. In particular, we can study classes of higher-order operators in this 
way as well. 



Apn(A) = a„p„+i(A) + /3„p„(A) + a„_ip„(A) 
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2 Preliminaries 

In this section we recall some results needed in the sequel. We first record the properties of the 
Laguerre polynomials. They satisfy the differential relations, [21, (4.6.13), (4.6.15)]: 



dx 



1 d 



x"e ^ dx 

A generating function of the Laguerre polynomials is 
f;LW(x)t" = (l-r"-'exp(^), 

n=0 ^ ^ 

see [21, (4.6.4)], [29, 37] and it implies 
The orthogonality relation is 



'x"e-^'L(;^)(x)4")(x)dx ^ r(a + n + 1) ^ 



n! 



m,n ) 



a > -1. 



The Meixner polynomials are, [21, § 6.1], [29, § 1.9], 
Mn{x; P,c) = 2Fi{-n,-x; f3;l - l/c), 
and have the generating function 



n=0 



V ^M„(x; (3, c)r = (1 - t/cTil - t)-^-^. 



The orthogonality relation is 



c " n! 



a;=0 



x! (/3)„(l-c)/5 



valid for /3>0,0<c<l and their three term recurrence relation is 

-xMn{x;P,c) = '-Mn+i{x;P,c) + M„_i(x;/3,c) 

1 — c 1 — c 



n + c(/3 + n) 
1 - c 



M„(x;/3,c). 



The Meixner polynomials satisfy the second order difference equation 



x\ ^/(/? + l),c- 



c^(/3)x. 



V 



x\ 



A Mn{x;(3,c) 



nc — I 



M„(x;/3,c), 



where 



(2.1) 
(2.2) 



(2.3) 



(2.4) 



(2.5) 



(2.6) 



(2.7) 



(A/)(x) = fix + 1) - fix), (V/)(x) = /(x) - fix - 1). 
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The Meixner-Pollaczek polynomials {Pn^\x; 0)} satisfy the orthogonality relation [29, (1.7.2)] 

I /_.(^*--)lr(A + i.)PF/«(.;*)pW(..0)d.= ^11=^^ 

for A > 0, < I?!) < TT, and the three term recurrence relation [29, (1.7.3)] 

(n + l)Pi+\(x; 0) + (n + 2A - l)Pi^\(a;; 0) = 2[x sin + (n + A) cos 4>]Pi^\x; </>), (2.8) 

with P^^\x;(f)) = 1, Pi^^^(x;(/)) = 2[j; sin + A cos 0]. 

We parametrize the independent variable x hy x = {z -\- l/z)/2 and given a function we set 
f{z) = f{x). The Askey-Wilson operator Vg and the averaging operator Ag are defined by, [21, 
§ 12.1], 

where e(x) = x = (z + l/z)/2. 

The Askey-Wilson operator is well-defined on H1/2, where 

:= {/ : f({z + l/z)/2) is analytic for < \z\ < q-"} . 

Let Tiw denote the weighted space L^( — 1, 1; w{x)dx) with inner product 

(/, 9)w := f{xUx)w{x) dx, ll/IU := (/, f)'J' 
and let T be defined by 

Tfix) ■.= -^^Vg{pVgf) (X), 

for f m Hi. We shall assume that p and t/; are positive on (—1, 1) and also satisfy 

(i) p{x)/Vl-x^ £ l/peL(-l,l), 

(ii) w{x) G L(-l, 1), 1/u; e L ( -1, 1; "^"^ 



:i-x^ 



The expression Tf is therefore defined for f £ Hi, and the operator T acts in "Hu,. Furthermore, 
the domain Hi of T is dense in since it contains all polynomials. The following theorem is 
due to Brown, Evans and Ismail [12]. 

Theorem 2.1. The operator T is symmetric inT-L^ and positive. 
The Al-Salam-Chihara polynomials are defined by [21, (15.1.6)] 



\tie'^ ,tie-'^ 

tlt2,0 



Pnix;ti,t2\q) = 3(p2 [ ^ ';^_n^ 1^1 

Their weight function is 

( (e^-^e-^-^g)oo/sing 
w{costl;ti,t2) := -, ^ ^ ^ ^ — ^ — , (2.9) 

and their orthogonality relation will be stated in (7.4). The generating function for the Al- 
Salam-Chihara polynomials is [21, (15.1.10)] 



n=0 



Sr (^1^2; g)n f t \ .\ {ttl,tt2;q)oc ^ 

UJ = (tie-, tie--- q).^- ^'-''^ 
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Theorem 2.2. Consider the three term recurrence relation in orthonormal form 

xpn{x) = a„+ip„+i(x) + bnPn{x) + a„p„-i(x), n > 0, a„ > 0, 6„ G M, (2.11) 

with aoP-i(a;) := 0. Then the moment problem is determinate, that is, it has a unique solution, 
if one of the following conditions hold 

5^^^ = oo, (2.12) 

an+lOn+2 

o-n + bn + On+i < C, for some C, (2.13) 
CLn — bn + On+i < C, for some C. (2.14) 

The condition (2.12) is Exercise 2 on p. 25 of [1], while (2.13), (2.14) are Theorem VII.1.4 
and its corollary in [11, pp. 505-506]. 

Theorem 2.3. Let pn{x) be generated by (2.11). Then the zeros of the polynomial pn{x) are in 
{A,B), where 

B = maxjxj : < j < n}, A = min{yj : < j < n}, 
where yj < xj and 



Theorem 2.3 is the special case Cn = 1/4 of a result due to Ismail and Li in [26]. The full 
result is also stated and proved in [21, Theorem 7.2.7]. 

The zeros of orthogonal polynomials are real and simple, so we shall follow the standard 
notation in [39] or [21] and arrange the zeros Xn,j, 1 < j < n as 

Xn,l ^ 3^)2,2 > ■ ■ ■ > Xrt^n- 

3 A differential operator related to the Laguerre polynomials 

Consider the differential operator 



iTLf){x) ^ ^ 



x°e ^ dx 



xdf 



dx 



(3.1) 



We will discuss a generalization of this operator in Section 6. The boundary value problem we 
are interested in is T^y = Xy with the boundary conditions x^^"^'^ f{x)e~^^'^ — >■ as x — >• and 
x — )■ cxD. The equation T^y = Xy is 

x'^y" + (a + 2)xy' - x'^y' = Xy. 

It is easy to see that Tl is symmetric on weighted L2 space with the inner product 

poo 

{f,g)= / x''e-^f{x)g{x)dx. (3.2) 
Jo 
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The [111,11) matrix elements of Tl as an operator in L2(0,oo,x'^e ^) in the basis {lI^\x)} 
can be calculated using (2.3), (2.4), (2.2); 



L 



(a+l) 
m— 1 



X)] 



dx 
d 



dx 



dx 



x°+2e-^^(4"+i)(.)-LtV^(x)) 



dx 



x^+'e'^[Li:+'\x) - L^^^l\x)][nLi-+'Hx) - (n - l)Lt\'\x)]dx 
r(n + a + 2) r(m + a + 2) r(n + a + 2) T{n + a + l) 



(n-1)! (m-1)! 
'5m+i,„ - (a + 2m) ^ 



r(m + a + 2) r(n + a + l) r(m + a + l) 



(n-1)! 



(n-2)! 



(m-iy. V- ' — ^ (m-1)! ' (m-2)! 

Thus the sought matrix representation of Tl is tridiagonal. It is also clear the constants are in 
the null space of Tf,, so we mod out by the constant functions. Let {Am,n{L)} be the matrix 

elements. Thus Am,niL) is 



(3.3) 



r(m+a+ijr(n+a+i) tim^s the above expression. Thus 

Am,n{L) = m^J (m + l)(m + a + l)(5rrt+l,n 

- m(2m + a)bm,n + (m - \)^Jm(a + m)(5m,n+i- 



The effect of modding out by the constants is to delete the first row and column of the matrix is 
to shift m and n by one. Thus we consider the tridiagonal matrix B = {Bm,n), m, n = 0, 1, . . . , 



B„ 



(m + l)A/(m + 2)(m + a + 2)6m,n-i 

- (m + l)(2m + a + 2)5m,n + m-y/ (m + l)(a + m + l)dm,n+i- 



Now the spectral equation EX = BX where X is a column vector, when written componentwise 
is a three term recurrence relations and the component of X are Pn{E) ■ The corresponding monic 
polynomials satisfy the three term recurrence relation 

Epm{E) = Pm+i{E) - (m + l)(2m + a + 2)pm{E) + n?(m + l)(m + a + \)p^-i{E). 

This is [29, (1.3.5)] and identifies the PmS as continuous dual Hahn polynomials with the 
parameter and variable identifications 

(a + 1)2 



1 — a 



1 + a 



3 + a 



E 



2 ' 2 ' 2 ' 4 

As stated in the Introduction the spectral decomposition is given in terms of the orthogonality 
measure for the continuous dual Hahn polynomials Sn{x;a,b,c), see [29, § 1.3]. The measure is 
absolutely continuous on [0, oo) in case a, b and c are positive. In case one of them is negative, 
finitely many discrete mass points have to be added. In the present case however we assumed 
a > —1, hence only a can be negative. Explicitly, with 5„(x) = Sn{x;a,b,c), 

2 



1 



27rT{a + b)T{a + c)T{b + c) 

M 



Sn{y )Sm{y ) 



r(a + iy)T{b + iy)T{c + iy) 



T{2iy) 



dy 



r(-2o)r(6 + c) ^ ^ ; ; ™v V ^ ^A;!(a)fc(a-6 + l)fc(a-c+l)fc^ ^ 



= 5nm.n\{a + b)n{a + c)n{b + c)n, 

where M = max{k £N : k + a < 0}. 

Now [24, Theorem 2.7] gives the following proposition. 



(3.4) 
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Proposition 3.1. The unbounded operator acting on the subspace of polynomials in 
L2(0, oo, x"e~^) is essentially self-adjoint. The spectrum of the closure has an absolutely conti- 
nuous part (— oo, — |(a + 1)^]. The discrete spectrum consists of {0} and {Ek \ k G {0,...,M}}, 
where M = max{k e N | A: + (1 - a)/2 < 0} and Ek = {k + l){k - a). 

The explicit spectral measure can be obtained from the orthogonality measure for the con- 
tinuous dual Hahn polynomials (3.4), cf. Section 1. 

This discussion of the differential operator Tl is related to the set-up of [23], where the case 
related to Jacobi polynomials is considered. In [23] we assume that we did not need to mod 
out a null space. The differential operator Tl can be related to the confluent hypergeometric 
differential equation in the same way as the hypergeometric differential operator shows up in 
[23, § 3]. 



4 A difference operator related to Meixner polynomials 

The generating function (2.5) implies 

/3Mnix;(3,c) = + n)Mnix; P + l,c) - nMn-iix; 13 + l,c). (4.1) 
The second order linear operator to be considered is Tm, 

(TMf){x) := -^V (^-^^^Af) (x). (4.2) 



c^(/3)x V a;! 

We consider the inner product spaces endowed with the inner product 

(/'5)/^=f;^/(^)^- (4-3) 

x=0 

The operator Tm is formally self-adjoint with respect to (•, 
Using (4.1) and (2.7) we see that 

TuMnix- P, c) = ^) W + n)Mn{x; /3 + 1, c) - n(n - l)M„_i(x; /? + 1, c)] . 

Therefore 

^'^^^(M^(x; /3, c),TMMn(x; /?, c))^ 
(c- 1) 

= ((/3 + m)M^(x;/3 + l,c) -mM„_i(x;/3 + l,c), 

n(/3 + n)Mn{x; /3 + 1, c) - n{n - l)Mn-i{x; (3 + 1, c))p+i 
= m{f3 + m(fhm{l3 + l)(5m,n - ^{m + l)(/3 + m)hm{l3 + l)5m,n~i 

- m{m - l)(/3 + m - l)hm-iiP + l)5m,n+i + m^(m - l)hm-ii/3 + l)dm,n, 

where 

C ^Tl^ 

hniP) = (M„(x; /3, c), Mn{x; /3, c))^ 



Since Tm annihilates constants we mod out by the space of constants and let the matrix elements 
of Tm be {Bm,n : m, n > 0}. Thus 

D {Mm+i{x](3,c),TMMn+i{x;[3,c))i3 

BmA^l) = ====== . 

hm+l{f3)hn+l{l3) 
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In other words 

c/3(/3 + l)Bm,n = -[{m + l){m + (3 + l) + m{m + l)c](5„,„ (4.4) 
+ my^ c(m + l)(m + I3)6m,n+i + im + l)v^c(m + 2)(/3 + m + l)(5m,„_i. 

Now scale the energy parameter Ehj E = —x/{cf3{f3 + 1)). This translates into the monic three 
term recurrence relation 

xPm{x) = P,n+i{x) + [(m + l)(m + /3 + 1) + m(m + l)c]P„(x) 

+ cm2(m + l)(/3 + m)Pm_i(x). (4.5) 

The polynomials generated by (4.5) seem to be new. They give a one parameter generalization of 
the continuous dual Hahn polynomials which is different from the Wilson polynomials. Finding 
the orthogonality measure of these polynomials remains a challenge. It clear that the measure 
is unique and is supported on an unbounded subset of [0, oo). They are birth and death process 
polynomials corresponding to birth rates 6„ = (n + l)(/3 + n + 1) and death rates dn = cn(n + 1), 
see (5.6). By (2.14) the corresponding moment problem is determinate. So in this case we do 
not have a precise analogue of Proposition 3.1, except that the unbounded operator Tm defined 
on the polynomials has a unique self-adjoint extension. 



5 Operators with additional potential 

We consider the case of second order operators which arise from classical orthogonal polynomials. 
Let Pn{x) be a monic family of classical orthogonal polynomials and T a second order operator 
such that 

Tpn{x) = XnPn{x). (5.1) 

Also assume that the three term recurrence relation for the pnS is 

XPn{x) = Pn+lix) + anPn{x) + /3„p„_i(x). (5.2) 

We now consider the spectral problem 

{T + -ix)il:{x,E) = Eiljix^E). (5.3) 

One can think of Ty = Ey as a free particle problem then (5.3) will be a Schrodinger problem 
with potential 7X. Let fi be the orthogonality measure of {pn{x)} and assume we deal with the 
case when the polynomials are dense in L2(M, //), which is true with very few exceptions. The 
orthonormal polynomials are {pn{x)/\/ l3iP2 ■ ■ ■ /3n } and form a basis for L2(M, /i). The matrix 
element of T + 7X with respect to this basis are 

1 1 /■ 

Bm,n = —^ n / Pm{x){T + 'yx)pn{x)dfl{x) . 



j=l k=l 



Clearly 



Bm,n = (An + 7"n)'5m,n + 7\//3m<5 m,n+l 

The matrix B = {-Bm,n} is tridiagonal and generates the monic orthogonal polynomials via 
{ME)] 

ME) = l, ME) = E - Xo - jao, 

(Pn+iiE) = {E-\n- 7an)ME) - j^Mn-iiE). (5.4) 



Spectral Analysis of Certain Schrodinger Operators 



9 



We can still scale E hy E = ^[x — rj) and introduce additional parameters to help identify the 
polynomials as known ones. Thus we let tpnix) = S,~"(f)n{E) and transform (5.4) to 

^o(a;) = 1, iJi{x) = X -Tj - (Xo + ^/ao)/^, 

Ipn+lix) = [x-r]- {Xn + 'yan)/C]'ipnix) - {-f / f3n1pn-lix) . (5.5) 

The importance of this scaling will be made clear in the examples. 

Recall that (5.2) generates a birth and death process polynomials if there are sequences {bn} 
and {dn} such that 

an = bn + dn, = C^n^n-l, (5.6) 

and for n > 0, &n-i > and dn > 0, with do ^ 0. One can represent the transition probability 
of going from state m to state n in time t as the Laplace transform of the product of two 
orthogonal polynomials and their orthogonality measure. For details, and additional information 
and references see [21, Chapter 5] and the survey article [25]. 

The hnS and dnS are birth and death rates at state (population) n. In the case of birth and 
death processes with absorption (killing) Karlin and Tavare [28] showed that the corresponding 
orthogonal polynomials satisfy (5.2) where (5.6) is now replaced by 

an = hn + dn + Cn-, /?„ = dn^n-l, 

where c„ is the absorption rate at state n. This leads to the following remark. 

Remark 5.1. Assume that T is a positive linear operator and (5.1) holds where {pn} are birth 
and death process polynomials with birth and death rates {bn} and {dn}, respectively. Then 
the orthogonal polynomials in (5.5) with ^ = 7 which arise in the tridiagonalization of T + 7X 
are polynomials associated with a birth and death process with absorption where the birth and 
death rates and {dn}, respectively and the absorption rates are {Xn/l}- 

The phenomena described in Remark 5.1 seem to be related to shape invariance and related 
topics in Discrete Quantum Mechanics recently developed by R. Sasaki and his coauthors, see 
the recent survey [35]. 

Example 5.2 (Laguerre polynomials). In this case T is as on the left-hand side of (2.2) and 

A„ = -n, a„ = 2n + a + l, /3„ = n(n + a), pn{x) = {-l)'^n\L^^\x). 
The recursion in (5.5) is 



V'n+Kx) 



X — 77 + 



71(1-27) 7(« + l)l,,. .^/.^ , „n7^ 



^n{x) - n{n + a)-^ipn-i{x), (5.7) 



with E = ^{x — rf). When 7 = 1/4 we take = 1/4, r/ = —2a — 2. This identifies ipn{x) 
as {—1)'^ L!>n\x) . Hence the spectrum is purely continuous and is given by x > 0, that is 
E > {a + l)/2. The absolutely continuous component is given by 



AE) = 4°+i ^"P (^(°- + l)) [j; _ + l)/2]-e-^^, Ee[{a + l)/2,oo). 
r(a + 1) 

We next assume 7 > 1/4 and compare (5.7) with the following monic form of (2.^ 

, , , , ,1 , / x n(n + 2A — 1) , , . 

Wn+i{x) = X - (n + A) coi<p\ipn[x) r— o- yJn-iKx). 

4 sm d) 
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We make the parameter identification 

7 = ^sec2(,/,/2), e = -tan((/>/2), A = (a + l)/2, = cot((/>/2). (5.8) 

With this choice of parameters we identify the ijj^s as Meixner-Pohaczek polynomials. Indeed 
'>Pn{x) = Pi'^\x; (p) where 

7 = (l + ^2)/4, A = (a + l)/2, ^ = -^, (5.9) 

and 

ME)=r''Pi''Hv + E/0- (5.10) 

The spectral measure /x is absolutely continuous and when normalized to have a total mass 1, 
its Radon-Nikodym derivative is 

2^ V^'-M24>--)-)\ri^, + ia + l)/2)\\ 



1 + 7rer(a + l) 

We now consider the case < 7 < 1/4. We identify (5.7) with the monic form of (2.6), namely 



Vn+lix) 



n — c(/3 + n) 



cn{n + /3- 1) 

yn{x) 7- ^2 yn-l{x). 

(1 - cY 



1 -C 

This is done through the parameter identification 

\/c f 1 - \/c /3Vc 



:i + ^)2' ' 1 + ^' ^ ' ' ' l-^c 

It is clear from (2.6) that the ?/„'s are monic Meixner polynomials. 

Note that such a division also occurs in the spectral decomposition of suitable elements in 
the Lie algebra su(l, 1) in the discrete series representations, see [31, 33]. So one can ask for Lie 
algebraic interpretations along the lines of [34], see [36] for a related result. 

Example 5.3 (Ultraspherical polynomials). In this case [21, (4.5.8)] 
T = (1 - x^r^^y^j- ((1 - x^Y^'l'£j , A„ = -n(n + 2.). 

The coefficients in the monic form of the three term recurrence relation are, see [29, (1.8.18)], 

n(n + 2v — 1) 
4(n + ujyn + V — 1) 

Thus the recursion in (5.5) becomes 

i:n+i{x) = [x - ?? + r^n{n + 2v)\Mx) - ^^^^i ~ !| V^n-i(x). (5.11) 

4(n + z^j(n + 1/ — Ij 

We do not know the orthogonality measure of the polynomials in (5.11). The special case v = \ 
appeared earlier in the work of Alhaidari and Bahlouli [6, (3.8)] where they applied the J-matrix 
method to quantum model whose potential is an infinite potential well with sinusoidal bottom. 
The same case also appeared in the work [18] by Goh and Micchelli on certain aspects of the 
uncertainty principle. Determining the orthogonality measure of these polynomials will be very 
useful. 
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The parameter r] can be absorbed in x, hence we assume rj = 0. In the notation of (2.11) 
n(n + 21/) 9 ^'^n(n + 2u — 1) 



C ' 4e2(n + i/)(n + z^-l)' 

In the case ^ > 0, since bn < and bn — ^n— i ^ 0; Xheorem 2.3 imphes that the smallest zero 
of Pn{x) is approximately + bn-i) — ^\bn — bn-i\- Hence 

Xn,n = -r^n{n + 2u) + 0{l). 

On the other hand is monotone decreasing if u > lor— l/2<i/<0 and monotone increasing 
if < u < 1. Using 

libn + bn^i) ± ^Vibn-bn-i)^ + 16al 

< ^{bn + bn-l) + ^\bn - 6n-l| + 2a„ < 2 max{ai , Coo} , 

where Ooo = lini an- Therefore 

n— >oo 

Xn,i < 2max{ai,aoo}- 

Thus the spectrum is unbounded below and is contained in (—00, 2 maxjai, ffloo})- It is important 
to note that pi{0) = 0, hence the right end point of the spectrum, being lim Xn 1 is positive. 

The case = 1 is the case when our starting point is the Chebyshev polynomials of the second 

3 I /„2 I 9 



kind. In this case a„ = oi for all n and the largest zero of P2{x) is ~| + y Q^f + | > 0. 

Example 5.4 (g-Ultraspherical polynomials). The weight function is supported on [—1, 1] and 
is given by [21, § 13.2] 

w{x; I3)dx = (^g2»e|^e-2»-e.gJ^ ^^' x = cos ^, /3 < 1. 
The second order operator is [21, § 13.2] 

In this case 

(l-g'^)(l-/?2g'^-i) 



an = 0, /3„ 



4(l-/3g")(l-/3g"-i)' 
and the recurrence relation in (5.5) gives 



X-V+ " \, (i-g")(i-/3V) 



tpn{x) 



7^(l-g")(l-/3^g"-^) 

4^2(1 _/3gn)(i_^^n-l) 



It is clear r] can be absorbed in x so we may assume rj = 0. We do not know any explicit 
formulas for the above polynomials. It is clear that they are orthogonal on an unbounded set and 
that Condition (2.12) is satisfied, hence the orthogonality measure is unique. As in Example 5.3 
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we can show the the spectrum is bounded above and unbounded below and estimate the largest 
and smallest zeros of pn (x) . In the present case 

(l-g)2^ ^ 4^2(1 _^^n)(l_^gn-l)- 

Here again the h^s are negative and decreasing in n for ^ > 0. A simple calculation shows 
that a„ increases with n if < /3 < (7 and decreases with n \i q < (3 < 1. Thus 



A := max{a„ : n = 1, 2, . . . } 



if < /3 < 
ai if g < /3 < 1. 



Therefore Theorem 2.3 shows that the smallest zero Xn^n satisfies 

Xn,n > ^{bn + &n-l) + ^{K " " 2a„ = 6„ - 2a„ > bn - 2A. 

Indeed Xn,n = bn + 0{1). To determine the other end of the spectrum note that 

V {bn - bn-iy + 16a2 <\bn- bn-l \ + 'iUn = 6n-l " &n + 4a„. 

Thus 

libn + bn-i) + ^V{bn-bn-i)^ + l6al 

< ^{bn + bn-i) + ^{bn-1 " &n) + 2a„ = 5„_i + 2a„ < 2yl, n > 0. 

Consequently the largest zeros Xn^i is < 2^. Therefore the spectrum of T + 72; is unbounded 
below and is contained in (— oo,2A]. 

Example 5.5 (Chebyshev polynomials). The Chebyshev polynomials of the first and second 
kinds are special ultraspherical polynomials and special q-ultraspherical polynomials as well. We 
will only discuss the polynomials {{7„(x)} but the reader can easily write down the corresponding 
formulas for the polynomials {T„(x)}. 

The UnS correspond to 1/ = 1 of (5.11) and the case /3 = g of (5.12). Thus we are led to the 
following systems of orthogonal polynomials 

rn+i{x) = [x -ri + ^ ^n{n + 2)]r„(x) - -^rn-iix), 



Sn+l{x) 



l-nc-l 



7' 



Here again we do not know any explicit representations or orthogonality measures for the polyno- 
mials {rn{x)} or {sn{x)}. Again condition (2.12) is satisfied for {r„(j;)} and {sn(x)}. Therefore 
the orthogonality measures of both families of polynomials unique. 



6 Adding a linear potential 

In this section we yet have a variation on the problems of potential introduced at the beginning 
of Section 5. We start with (5.1) where the eigenfunctions satisfy (5.2). We then consider the 
Schrodinger operator 



S = {x + c)T + 7x. 



(6.1) 
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We illustrate this idea by considering the operators Tl and Tm for the Laguerre and Meixner 
polynomials defined in (3.1) and (4.2), respectively. 
The Laguerre case. Here we take c = and 



T 



d 



d^ 
dx 



x"+i dx 

With the notation in (3.1) we have 
S = Tl + ^x. 

We use the inner product (3.2) and our weighted L2 space is L2(0,oo,x"e~ 
elements Sm n are 



The matrix 



mini 



r(a + m + l)r(a + n + l 

I ^ 
dx 



X 



Li^Hx)^lx''+h''' 







d 1 f°° 

— 4")(x) Ux + 7 / L(^)(x)x4")(x)x"e-^dx 
ax J Jo 



Using the recurrence relation 

x4")(x) = -(n + l)4"i(x) - (n + a)Li%ix) + (2n + a + 1)4" )(x), 

and the calculation of the matrix elements in (3.3) we find that 

Sm,n = [7(2"i. + a + 1) - m{a + 2m)]5m,n 

+ (m - 7) v^(m + l)(m + a + l)5m,n-i + (m - 1 - -f)y/m{m + a)5m,„+i. 

The null space of 5" is trivial so there is no need to mod out by the null space as we did in 
Sections 3 and 4. The monic polynomials {pn{E)} which arise through tridiagonalization are 
generated by po{E) = 1, pi{E) = E — 7(0 + 1), and 



Epn{E) = Pn+i{E) + n{n + a){n-l- 7) V-i(^) 
+ [7(2n + a + l)-n{a + 2n)]pn{E). 



(6.2) 



The polynomials in (6.2) form a two parameter subfamily of the continuous dual Hahn polyno- 
mials [29, § 1.3] with the parameters 



-7-(a + l)/2, 



c = (a + l)/2, E 



1 



(a + 1)^ 



In the above analysis we assumed a > — 1, hence 6 = c > 0. If 7 < — (a + 1)/2 then S has purely 
a continuous spectrum supported on (—00, —{a + 1)^/4], see (3.4). If 7 > —{a + l)/2 there is 
also discrete spectrum, and we obtain as in Section 3 the following proposition. 

Proposition 6.1. The unbounded operator S acting on the subspace of polynomials in 
L2(0, 00, x"e~^) is essentially self-adjoint. The spectrum of the closure has an absolutely con- 
tinuous part (—00, — ^(a + 1)^]- The discrete spectrum consists of {Ej. \ k S {0, . . . , M}}, where 
M = max{A; G N | k--f-{l+a)/2 < 0} and Ek = {k--f){k--f-a-l) G (-|(a+l)^ 7(7+0+1)] . 

Remark 6.2. The case 7 = of Proposition 6.1 reduces to Proposition 3.1 using a similar 
reduction as in [23, Remark 3.4]. Note that in this case Eq = 0. 
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Remark 6.3. In this case S = T^ + jx can be written as 5 = x (1)°"'"^ + 7) , where is the 

second order differential operator 



D 



a+l 



1 



d 



for which L>"+iLi"+^^(x) : 



d 



dx 

nL^n~^^\x). This shows that S is of the type considered in [23]. 
The Meixner case. With the notation in (4.2) we let 



This corresponds to c = in (6.1), and T^j as in (4.2). We use the inner product (4.3) and 
our space is now L2 weighted with the orthogonality measure of the Meixner polynomials with 
parameters c and f3. The matrix elements Sm,n are 



(/3)^(/3)„(l-c)2/3 
(l-c)7 



^ M„(x; (3, c)V ( (^±^AM„(x; /3, c] 



x=0 



+ 



c/3(/3 + l) 



n E ^-^^rn{x; /3, c)xM„(x; /3, c) 



x=0 



We already calculated the matrix elements of Tm in (4.4), but we must replace m, n by m — 1, 
n — 1, respectively. Using the recurrence relation (2.6) we then compute the matrix elements of 
a constant times x. This leads to 

c/3(/3 + l)Sm,n = -[m{m + /?) + m(m - l)c + 7m + c7(/3 + m)]5r„,„ 

+ (m - 1 - 7) 1/ cm(m + /3 - 1) 6m,n+i + {m - 7)v^c(m + l)(/3 + m)5m,n-i- 

Therefore the corresponding orthonormal polynomials Pn{E) satisfy the three term recurrence 
relation 

c/3(/3 + l)Epm{E) = -[m{m + /3) + m(m - l)c + 7m + c7(/3 + m)]pm{E) 



+ (m - 7) c(m + 1) (/3 + in)pm+i (E) + (m - 1 - ^)^J cm(ra + /3 - l)pm-i(-E')- (6.3) 

We assume 0<c<l,/3>0, since we are dealing with the Meixner polynomials. In order 
to have (6.3) satisfy the conditions for an orthonormal polynomial system we assume 7 < 0. 
The polynomials generated by (6.3) seem to be new. Here again we have neither explicit rep- 
resentations or generating functions, nor do we know their orthogonality measure. One can say 
however that their orthogonality measure is unique since condition (2.13) is clearly satisfied 
for sufficiently large n. Using Theorem 2.3 and some estimates we see that the support of the 
orthogonality measure is contained in (— cxd, a] for some a > 0. 



7 The Al-Salam— Chihara polynomials 

Recall that w is defined in (2.9). The generating function (2.10) implies 

(1 - tit2)pn{x;ti,t2) = (1 - ht2q"')pnix;ti,qt2) - ^1*2(1 - q"')pn-i{x;ti,qt2) 
We will first consider the case when the operator T is 



(7.1) 



L :-- 



w{x;ti,t2) 



w[x 



■,q'/Huq'/h2)V, 
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Apply (15.1.6) and (12.2.2) in [21] to see that 

2tigi-"(l -g^ 



T>qPn{x;ti,t2) 



{l-q){l-ht2) 



(7.2) 



Moreover 



Vq[w{x,^/qtl,^t2)pn-l{x]^/qtl,^/qt2)] = ^^^1*"^^ W{x,tl,t2)pn{x]tl,t2). (7.3) 

For real ti, t2 with |ti|,|i2| < 1 the Al-Salam-Chihara polynomials {pn{x]ti,t2)} satisfy the 
orthogonality relation [21, (15.1.5)], 

hnitl,t2)Sm,n = j Pmix] h, t2)Pnix; h, t2)w{x; ti, t2)dx , (7.4) 

" ^' ^ {q,tit2;q)oo{tit2;q)n' 

and are complete in L^(— 1, l;w{x; ti,t2)dx). In view of (7.4) the orthonormal Al-Salam-Chihara 
polynomials are 



Pn{x;ti,t2) 



(9,^1*2; 9)00(^1*2; g)n 



Pn{x]ti,t2). 



Lemma 7.1. Let {Am,n{AC)} he the matrix elements of L in the basis {pn{x',ti,t2)} ■ Then 
4gi-™(i _ , 



+ 



+ 



4t2g^~™(l - q"" 
4t2g^"'"(l -9™"^ 



(1-9) 

Proof. Clearly equation (7.2) implies 

Pm{x; ti,t2)w{x; ti,t2)Lpn{x; ti,t2)dx 



[l-M2g™ + ti(9-g™)]<5„,„ 

V(l - g™+l)(l - tit2q"')6m,n-l 
2 V(l - g™)(l - tit2q"'-^)6m,n+l- 



2ti{l-q^)q 



Pm{x; ti,t2)Vq w{x] q^^'^ti, q^/'^t2)pn-i{x] tiy/q, t2^/q) 



(l-g)(l-iit2) 7-1 
In view of (7.1) the integrand in the last step is 



(1-M2<7™) I , tMl-q^ I , 

Pm{x; ti,qt2) ^ — Pm-i{x; ti, 9*2) 



dx. 



1 - tlt2 



w 



1 - tlt2 
(1 - tlt29") 



(x; q'/h,,q'/h2) { ^\ ^ Pn-l (x; v^ti, 9^/^*2) 



gti*2 

p„_2(x; y/qti,q^^'^t2) 



1 - gtit2 

Applying (7.3) we see that the quantity after the x is 
2 



hiq-l) 



w{x; ti,qt2) [(1 - tit2q"')pn{x; ti, 5*2) - ^1*2(9 - q^)pn-i{x; ti, 9*2)] • 
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Therefore 

TT-, „\ 1^ n / Pm{x;ti,t2)w{x;ti,t2)Lpn{x;tl,t2)dx 

w{x,ti,qt2)[{l - tit2q"')pm{x;ti,qt2) - tit2(l - q"^)Pm-i{x]ti,qt2)] 

-1 

X [(1 - tit2q"')pnix;ti,qt2) - tit2{q - q'^)pn-iix;ti,qt2)]dx 
= [(1 - ht2qn^hm{h,qt2) + 44(1 - q^){q - q^hrn-lih, qt2)]Sm,n 

- qtit2{l - g")(l - ht2q"')hra{h,qt2)6^,n-l 

- tMl - q"'){l - tit2q"'~^)hm-l(tl,qt2)6m,n+l- 

The result follows since 

, t2)w{x; ti,t2)Lpn{x; ti,t2)dx. 



1 /•! 

■Am,n — 1 11 \T / \ I Pm{x',ti,'t 

V lT'm(tl,t2)nn[tl,t2) J -1 



The monic orthogonal polynomials generated by the matrix jj^^A, with A = {Am,n{AC)} 
satisfy 

Pn+iix) = [x- q-^il - g"+i)[l - tif2(?"+' + 4q{l - q'')]Pn{x) 

- tlq^-^^{l - (?"+i)(l - g")2(l - ht2qnPn-i{x). (7.5) 

Note that the recurrence relation (7.5) is invariant under g — )• l/g after scaling and renaming 
the parameters. The recurrence coefficients grow exponentially, and by [11, Theorem VII. 1.5] 
we can easily check that the moment problem corresponding problem does not have a unique 
solution (indeterminate). Hence L with domain the polynomials is not essentially self-adjoint. 
Nothing is known about the explicit formulas of the polynomials generated by (7.5) or any of 
their orthogonality measures. 

The ^"■'■-Hermite polynomials. We now study the g^^-Hermite polynomials of Askey [9], 
Ismail and Masson [27]. They are generated by ho{x\q), hi{x\q) = 2x, and 

hn+i{x\q) = 2xhn{x\q) - g~"(l - g")/i„_i (xjg). 

Here we use the parametrization x = sinh^. Recall the definitions [21, Chapter 21], [22] 

/V/2e«) - /V'/'e«) 



/(x) = /(e«), (VJ):-- 



(gV2 _g-i/2)(e5 + e-«)/2' 



{AJ){x) = l[f{qy'e^)+f{q-^/'e% 

The corresponding moment problem is indeterminate but all the A^-extremal measures have 
been determined in [27], see also [13, § 4] for another proof. They are purely discrete and are 
enumerated by a parameter a G {q, 1). The support is {xnia) : n = 0, ±1, ±2, . . . } and 

/^ Ir-n/ ni / MA a^"^"^^™-!) ^ ^2 2n) 

Xn{a) = -[q "/a - ag"], K{xn{a)}) = 2 _ / 2 . ^ ' (7-6) 

z [ a , q/ a , q, q)oo 

where /U is the corresponding normalized orthogonality measure. The orthogonality relation is 
[ hm{x\q)K{x\q) dii{x) = q)J^^^. 
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The lowering operator is 
The second order operator equation satisfied by the g~^-Hermite polynomials is [22] 

gV2(l + 2x2)p2y + Ji_^^^p^y = XnV, Xn ■= - ^ ^ " (7-7) 

With the measure fi defined as in (7.6) the matrix elements of the operator Th + 7X with Th 
the operator on the left side of (7.7) on L2{M.,fj,) with basis = q"'^"'~^^^/^hn{x\q) / \J {q; q)n} 
are given by 

{Th + ix)hn = |g-("+l)/Vl-(?"+l/ln+l + KK + l^q-""'^ ^/l^K-L 

Therefore the corresponding monic polynomials {p„(-E)} are generated by 

4g(l-g")^ 72 
_ 

This is essentially a perturbation of the Jacobi matrix of the g~^-Hermite polynomials by the 
diagonal matrix — ^^[^^^2 ^ • Apart from the shift ^j^^j^Id, this is a compact perturbation of the 
Jacobi matrix for the g~^-Hermite polynomials. Using [11, Chapter VII, Theorem 1.5] the mo- 
ment problem corresponding to the orthogonal polynomials generated by (7.8) is indeterminate. 
We now give bounds for the zeros oi pn{E). In the present case 

.2 l\.-nn „n^ , 4g(l-g") 



4^ ^ ^ ^' (l-g)2 • 

Take 7 > 0, it is clear that an is monotonic increasing while bn is monotic decreasing. The use 
of 

V {bn - 6n-i)2 + I6al < \bn - + 4a„, 
shows that the Xj's and y^-'s in (2.15) satisfy 

< + bj^i) + -bj) + 2a j = bj-i + 2a j < 2an, 

Vj > 7;ibj + bj-i) - l.{bj-i - bj) - 2aj = bj - 2aj > -2a^ 



2V"J — J-^/ 2'^"" ' ' " (l-g)2' 

Therefore 

4g 



Xn,l < 7g""/Vl-9". Xn,n > -79""^^^/^" 
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